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In recent years, searches of archival X-ray data have revealed galaxies 
exhibiting nuclear quasi-periodic eruptions with periods of several hours. 
These are reminiscent of the tidal disruption of a star by a supermassive 
black hole. The repeated, partial stripping of a white dwarfin an eccentric 
orbit around an -10° M, black hole provides an attractive model. A separate 
class of periodic nuclear transients, with much longer timescales, have 
recently been discovered optically and may arise from the partial stripping 
of amain-sequence star by an -107 M, black hole. No clear connection 
between these classes has been made. We present the discovery of an X-ray 


nuclear transient that shows quasi-periodic outbursts with a period of 
weeks. We discuss possible origins for the emission and propose that this 
system bridges the two existing classes outlined above. This discovery was 
made possible by the rapid identification, dissemination and follow-up 
ofan X-ray transient found by the new live Swift-XRT transient detector, 
demonstrating the importance of low-latency, sensitive searches for 


X-ray transients. 


Swift JO023017.0+283603 (hereafter Swift JO230) was discovered in Swift 
X-ray Telescope (XRT) data by the Living Swift-XRT Point Source (LSXPS) 
catalogue’s real-time transient detector! on 2022 June 22 (ref. 2). The 
source was Serendipitously present in an observation of an uncon- 
nected source, SN 2021afk (4.3’ away) and had a 0.3-10 keV count rate 
of 2.706 x 10-? counts s”. This field had been observed on 11 previ- 
ous occasions by Swift between December 2021 and January 2022. 
After combining all those observations, Swift J0230 was still undetected 


downto a3oupper limit of 1.5 x 10° counts s”. The last of these obser- 
vations was 164 days before the discovery of Swift JO230, placing a 
rather loose lower limit on its switch-on time. For convenience, we give 
alltimes relative to modified Julian date (MJD) 59752 (midnight onthe 
day of the discovery). The best localization of Swift JO230 is from the 
XRT, being RA = 02 h 30 min 17.12 s, dec. = +28° 36’ 04.4” (J2000), 
with an uncertainty of 2.8” (radius, 90% confidence). This is consistent 
with the nucleus of the galaxy 2MASX JO2301709+28360S0, but also 
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Swift JO230 


Fig. 1| Location of the new transient, Swift J0230, relative to its host galaxy 
and an old supernova. The image is an archival Pan-STARRS“ image of 2MASX 
JO2301709+2836050, with colour scaled arbitrarily for aesthetic purposes. The 
broken circle shows the 90% confidence Swift-XRT position of Swift J0230, and 
the solid one SN 2020rht. 


marginally consistent with the type-II supernova SN 2020rht (3.1” away), 
which had been discovered 2 years earlier on 2020 August 12 (ref. 3; 
Fig. 1). An optical spectrum of 2MASX JO2301709+2836050, obtained 
with the Nordic Optical Telescope (NOT; Fig. 2), gives a redshift 
z= 0.03657 + 0.00002. Assuming standard cosmological parameters’, 
this corresponds to a luminosity distance D, = 160.7 Mpc. The galaxy 
type is unclear, but it is either quiescent or, at most, a very weak active 
galactic nucleus (AGN; Supplementary Information). 


Results 

X-ray analysis 

Following the initial, serendipitous discovery, we undertook regular 
monitoring with Swift (see Methods for details). The initial outburst 
continued for 4 d following the discovery. On the fourth day, it ended 
with a rapid decay, the luminosity falling by a factor of 20 in just 57 ks. 
There was a brief rebrightening (by a factor of 4.5 in 6 ks), before it 
became too faint to detect. The light curve is shown in Fig. 3. Fitting 
this decay of luminosity with a power law, L œ (t — to)“ (where ty is the 
start of the first bin in this observation), gives a = 11.0 + 1.7. Eight sub- 
sequent outbursts were observed at ~25 d intervals, with durations of 
~10-15 d. The fifth outburst was either substantially longer 
(up to ~32 d), or consisted of a weak outburst, a return to quiescence, 
and then a second, longer outburst. This was followed by along gap of 
~70 d during which two possible short and weak outbursts were seen, 
before another outburst like the early ones. 

A Lomb-Scargle analysis (Methods and Extended Data Fig. 1) 
revealed moderately significant peaks at approximately 22 and 25 d 
periods, although each peak is ~1 d wide, confirming the quasi- rather 
than strictly periodic nature of the variability, as may be expected 
from eyeballing the light curve. Further consideration of the vari- 
ability required us to define what constitutes an outburst. The initial 
outburst and that from days ~41-48 appear clearly defined, but during 
the outburst from days ~60-75, the source underwent a sudden decline, 
being undetected on day 72 with an upper limit of L < 8.9 x 10" erg s” 
(0.3-2 keV), recovering tol ~ 2 x 10” erg s ‘by day 74. It seems plausible 
to interpret this as a single outburst with a sudden, brief dip. Hereafter, 
what constitutes an outburst becomes more subjective. The outbursts 
starting on days 89 and 102 could each be explained as comprising two 
short outbursts close together or a single outburst with a quiet phase 
in the middle. Note that in the first of these, if we sum the upper limits 
during this quiet phase, we find a detection at a higher level than the 
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Fig. 2 | Optical spectrum of the host galaxy 2MASX J02301709+2836050 
obtained with the NOT on day 132. Top, the black line shows the observed 
spectrum, while the orange line shows the fit to the stellar continuum provided by 
STARLIGHT. The vertical lines mark prominent emission and absorption features, 
which together allowed us to measure the redshift z = 0.03657 + 0.00002 
(Loconfidence). Bottom, the residuals between the observed data (stellar + 
nebular spectrum) and the fit (stellar continuum), which single out the nebular 
emission. The emission-line fluxes were measured from the residual spectrum 
and allowed us to place the galaxy ona BPT plot (Extended Data Fig. 3). 


upper limits found in the quiescent phases. During the long, largely 
quiescent period from days 111 to 195, the source was twice briefly 
detected with L ~ 1-2 x 10*' erg s7, but these are hardly ‘outbursts’ in 
the same way as the earlier emission. Based on visual inspection of 
the light curve, we defined an outburst as comprising any times where 
L>2x10“ erg s”. Details of the outbursts thus identified are given 
in Table 1. See Methods for full details of how these were derived. In 
summary, we detected transient X-ray emission that rapidly switches 
onand off again with a recurrence timescale that is of the order of 25 d 
but whichcan vary by several days between outbursts. The duration of 
the outbursts also shows significant variability with the longest being 
of the order of 20 d and the shortest less than a day. 

The X-ray spectrum during the outbursts was very soft, with no 
emission seen above 2 keV, and could be modelled well with a simple 
black-body emitter with only Galactic absorption. Due to this soft spec- 
trum, the typical energy bands used for XRT hardness ratios were inap- 
propriate. We selected 0.3-0.9 keV and 0.9-2 keV as this gave roughly 
equal counts in the two bands, thus maximizing the signal-to-noise 
ratio. The time evolution of this hardness ratio (Fig. 3) shows a clear 
correlation between luminosity and spectral hardness (Spearman 
rank P value of 1.3 x 10 for the data being uncorrelated), ruling out 
a change in absorption as the cause of the flux variation. We fitted 
the absorbed black-body model to each observation in which Swift 
JO230 was detected. The black-body temperature obtained is strongly 
correlated with luminosity (P value, 4.5 x 10%; Extended Data Fig. 2). 
No evidence for absorption beyond the Galactic column was seen. 

As noted earlier, although coincident with a galaxy nucleus, the 
XRT position is also potentially in agreement with that of SN 2020rht. 
We obtained a3 ks Chandra director’s discretionary time observation 
during the fourth outburst to obtain a better position, but unfortu- 
nately, this observation fell on day 97, which turned out to be in one of 
the mid-outburst quiet phases. 


Optical and ultraviolet analysis 

At optical and ultraviolet wavelengths, there is no evidence for out- 
bursting behaviour. We obtained data from both the Swift Ultraviolet/ 
Optical Telescope (UVOT) and the Liverpool Telescope (LT; Extended 
Data Tables 1 and 2). The host galaxy, 2MASX JO2301709+2836050, is 
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Fig. 3 | The temporal evolution of Swift J0230. All error bars are lo significance. 
Left, X-ray time series, binned to one bin per observation. Top, 0.3-2 keV 
luminosity light curve. The red data points marked with dots are the Swift-XRT 
detections. Grey arrows are 30 upper limits from XRT. The dark blue upper limit 
marked with an asterisk is from Chandra. The broad bins marked with crosses 
were created by combining consecutive XRT non-detections (upper limits in cyan 
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and detections in green). See Methods for full details. Bottom, the (0.9-2)/ 
(0.3-0.9) keV hardness ratio. The spectral hardness is strongly correlated with 
the luminosity. The vertical bands are at 25 d intervals. Right, light curve of 
the XRT observations taken on days 4 and 5, with one bin per spacecraft orbit, 
showing the rapid decay at the end of the first outburst. All error bars are at 
the lolevel. 


Table 1| Constraints on the timescales of the observed outbursts 


Outburst Start End Duration Period 
MJD min MJD max MJD min MJD max min (d) max(d) min (d) max(d) 

1 59,587.9 59,752.3 59,757.3 59,757.4 4.9 169.5 

2 597763.2 59,773.0 59,7731 59,787.5 01 24.3 10.9 185.1 
3 59,787.5 59,793.7 59,799.8 59,800.4 6.2 12.9 14.5 30.4 
4 59,811.0 59,815.3 59,826.4 59,827.0 1.1 16.0 WES) 278 
5 59,839.2 59,841.5 59,862.3 59,867.5 20.8 28.3 24.0 30.5 
6 59,870.8 59,886.7 59,886.7 59,887.9 0.0 17.0 29.3 47.5 
7 59,908.3 59,909.7 59,909.7 59,910.8 0.0 2.4 21.6 38.9 
8 59,942.3 59,948.0 59,954.3 59,955.6 6.4 13.4 32.5 39.6 
9 59,960.2 59,963.0 59,971.7 59,974.8 8.7 14.7 12.2 20.7 
10 60,010.8 60,011.9 60,022.4° = 11.6° = 478 51.8 


“The final outburst was ongoing when observations ended. 


clearly detected in all observations, but there is no evidence for vari- 
ability or an increase in flux compared to catalogued values. We also 
analysed UVOT data from the prediscovery observations and found 
no secure evidence for a change in brightness between those data 
and the observations taken during the outburst. Full details are given 
in Methods. 


Discussion 

The peak luminosity of -4 x 10” erg s”, the timescales of the outbursts 
and their quasi-periodic, quasi-chaotic nature, the soft X-ray spectrum 
and the lack of optical variability place strong constraints on possible 
models that explain Swift JO230. Although the lack of detection with 
Chandra means that we cannot rule outa positional association with SN 
2020rht, it is difficult to see how a supernova could have evolved into 
the object that we have detected. The spectrum, luminosity and vari- 
ability timescale are inconsistent with the properties of ultra-luminous 


X-ray sources’. Although certain supernovae could in principle be fol- 
lowed by X-rays froma newly formed millisecond pulsar‘, this should 
occur while the supernova is still visible in the optical. Margalit et al.’ 
showed how this emission could be delayed, but the timescales and 
luminosities they predict (for example, their Fig. 5) are not consistent 
with our observation. Equally, neither model explains the variability 
or spectrum we see in Swift JO230. We discuss this further in the Sup- 
plementary Information. 

We suggest that (near) periodic mass supply into an accre- 
tion flow onto the central supermassive black hole in 2MASX 
JO2301709+2836050 is the most likely explanation for Swift JO230. 
From simple energetics (Supplementary Information), we can infer 
that the total mass accreted during a typical outburst is -10% M,. In 
an AGN, flares and outbursts are common froma supermassive black 
hole at the heart of the galaxy accreting froma surrounding disc of gas. 
However, the timescale and spectrum of Swift JO230 are not consistent 
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with typical AGN behaviour, and 2MASX JO2301709+283605S0 itself 
does not appear to be an AGN (Extended Data Fig. 3; see Supplementary 
Information for a full discussion). We, therefore, consider the possibil- 
ity that one (or several) stars are interacting with, and feeding mass 
onto, the central supermassive black hole. One possible mechanism for 
producing the mass flowis the interaction of two stars in orbit around 
the black hole. If they pass sufficiently close to each other, material 
can be liberated from one or both stars that can feed the central black 
hole®’. To generate the required timescales from this model requires a 
pair of stars orbiting in the same direction and in the same plane’. This 
could occur for stellar orbits that are initially randomly oriented, if they 
can be subsequently ground down into the same plane by interaction 
with an AGN disc”. For Swift 0230, which lacks any clear signature of 
a standard AGN disc, it is unlikely that any stars orbiting the central 
black hole have the required orbits to achieve the observed timescales. 

Another possibility is a repeating, partial tidal disruption event 
(rpTDE), in which a star ona bound, highly eccentric orbit loses some 
of its envelope every pericentre passage due to tides from the black 
hole’s gravitational field. These events are a subclass of TDEs in which 
the ‘regular’ scenario sees the incoming star approach the black hole 
ona parabolic orbit and the star is destroyed by the first encounter 
(see refs. 11,12 for reviews of TDEs). 

The rpTDE model was investigated before the discovery of any 
corresponding sources (for example, ref. 13) and was suggested as the 
explanation of X-ray flares in the active galaxy IC 3599 (ref. 14). Recently, 
it has been proposed (for example, refs. 15-17) as a possible explanation 
for hours-long quasi-periodic eruptions (QPEs) discovered in galactic 
nuclei (for example, refs. 18-22). These works focused on the possibility 
of awhite dwarfinteracting with a relatively low-mass central black hole 
of mass ~10* * M,. A second set of sources show much longer outbursts 
(both in duration and recurrence period) and have been referred to 
as periodic nuclear transients (PNTs). These may be the same rpTDE 
mechanism acting with a main-sequence star, rather than a compact 
star”, and a more massive black hole (-10” § M,). 

In the rpTDE model, the donor star is in a highly eccentric orbit 
around a black hole. At each pericentre passage, the star has to 
approach, but not quite reach, the tidal radius at which the star would 
be fully disrupted. The outer layers of the star are liberated, and some of 
this material accretes onto the central black hole to power the outburst. 
The recurrence time of the outbursts is related to the orbital period of 
the star. The majority of the energy released from the accretion process 
occurs in the central regions near the black hole where the matter is 
most likely in the form ofan accretion disc. We can, therefore, provide 
an estimate of the black hole mass in Swift JO230 by comparing the 
temperature of -100 eV (-10° K), measured from the X-ray spectrum, 
with the peak temperature of a standard accretion disc”. This yields a 
mass estimate for the black hole of -2 x 10° M, (Supplementary Informa- 
tion). This is close to the mass estimates for QPE sources (for example, 
refs. 18,19,22). Note that QPE sources and Swift JO230 show very little in 
the way of optical emission, whereas the PNTs show strong optical emis- 
sion. This may simply reflect the difference in the black hole masses. 
That is, the soft X-ray spectrum and the lack of optical emission seen for 
Swift J0230 appear consistent with this estimate of the black hole mass. 

If the accreted material is stripped from the orbiting star during 
the pericentre passage of a highly elliptical orbit and the pericentre 
distance is a few gravitational radii (required to liberate any material 
fromthe surface of a white dwarf for black hole masses of a few 10° M,), 
then we would expect the outburst duration to be similar in different 
systems, regardless of their orbital (and hence outburst) period. This 
is because the pericentre passage of a highly elliptical orbit is approxi- 
mately that of a parabolic orbit and its duration is not connected to 
the orbital period. This means that it is difficult to explain both Swift 
JO230 (outburst duration of days) and the pre-existing QPEs (outburst 
duration of hours) as anrpTDE of a white dwarf around a modest-mass 
black hole. On the other hand, an rpTDE of a main-sequence star by 


an ~10’-108 M, black hole has been proposed to explain the PNTs 
ASASSN-14ko, for which the recurrence timescale is 114 d (refs. 23,24), 
and AT2018fyk, which exhibited a significant rebrightening after 
around 600 d of quiescence”. Swift JO230 clearly lies between these 
two classes of object. 

An important question is how the star arrived on such an orbit 
around the central black hole. Tidal capture, in which an orbiting star 
loses orbital energy due to tidal forces and becomes bound to the black 
hole’®, is typically incapable of generating the required orbits; however, 
the Hills mechanism” has been proposed asa viable formation route for 
the PNT ASASSN-14ko (ref. 30). In this mechanism, a binary star system 
approaches the black hole witha small enough pericentre distance such 
that the tidal force from the black hole is stronger than the gravitational 
force holding the binary together. This results in the binary being 
disrupted, with one component being ejected from the system and 
the other becoming locked into a bound, but highly eccentric, orbit 
about the black hole. If the progenitor of Swift JO230 were a binary 
consisting of a low-mass, main-sequence star and, say, a white dwarf, 
then the main-sequence star needs to be captured into a bound orbit 
around the black hole with the observed ~25 d period. For a black hole 
mass of M. ~ 4 x 10° M, this period corresponds to the most likely out- 
come of Hills capture from sucha binary system (the calculations of *° 
show that higher black hole masses are allowed but are significantly 
less likely to result in this period for the bound star; see Supplemen- 
tary Information for details). This is consistent with the mass estimate 
derived from the temperature of the X-ray spectrum. Further, the 
expected accretion timescales from such a system (Supplementary 
Information) are also consistent with those observed in Swift JO230. 

The variable shape and timescales of the outbursts seen in Swift 
JO230 may also be explained by this model. In a standard TDE, as 
opposed to anrpTDE, the star arrives on a parabolic orbit, meaning 
that some of the stellar matter is bound to the black hole (the inner tidal 
stream), and this material forms an accretion flow, whereas the rest of 
the stellar debris (the outer tidal stream) is unbound and leaves the sys- 
tem”. InanrpTDE, the star must be ona bound orbit around the black 
hole. In this case, both the inner and outer tidal streams can remain 
bound to the black hole. The inner stream falls back soonest and, thus, 
has a higher mass return rate, whereas the outer stream can return on 
longer timescales. Due to relativistic precession of the stellar debris 
orbits (both apsidal and nodal), the returning streams can collide and 
partially cancel their angular momenta to augment the accretion rate 
onto the black hole, with the magnitude of the effect depending onthe 
orientation of the colliding orbits (see refs. 31,32 for similar variability 
induced in accretion discs due to relativistic precession; the processes 
described therein may also occur in the discs formed in Swift JO230). 
The exact details of this interaction between the two streams, the accre- 
tion flow and the orbiting star are complex and require a full numerical 
analysis, whichis beyond the scope of this discovery paper; however it is 
clear that such an interaction will produce variable emission that could 
at least partially erase the more exactly periodic nature of the stellar 
orbit. An example of such complex interacting debris streams can be 
seen in Fig. 9 of ref. 33. It will be particularly important to determine 
if the star itself can be sufficiently perturbed during each pericentre 
passage, with, for example, tides imparting variations in oscillation 
amplitudes and rotation frequency and, thus, change the amount of 
mass transferred and the structure of each outburst. Additionally, the 
sharp decline observed at the end of each outburst may be driven by 
the returning star disturbing the accretion flow. These questions can 
be addressed with future theoretical investigations. 

We have proposed a single explanation for the QPEs and PNTs, 
as a repeated, partial, tidal disruption of a star in an eccentric orbit 
around a supermassive black hole. We have reported the discovery of 
an object that can bridge the gap between these classes. The QPEs are 
thought to be due to a white dwarf and a modest black hole (-10° M,), 
whereas the PNTs are thought to occur for a main-sequence star and 
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amore massive black hole (~10’ M,). Swift JO230 represents an inter- 
mediate class of system, which is consistent with a main-sequence 
star orbiting amodest-mass black hole. Given the timescales, modest 
fluxes and lack of emission outside of the X-ray band, Swift JO230-like 
systems are difficult to discover. Unlike QPEs, which can be discovered 
in archival data, their timescales and behaviour are not exposed by a 
single observation. It is only with the recent creation of a real-time 
transient detector’ that objects like this can be found rapidly enough 
for follow-up observations to expose their behaviour. That this event 
was found within 3 months of enabling this real-time search suggests 
that they are reasonably common, and we can expect to discover more 
objects of this class with sensitive, wide-field X-ray instruments such 
as eRosita™ and in the near future, the Einstein Probe”. 


Methods 

Discovery 

At 13.58 UTC on 2022 June 22, the LSXPS real-time transient detector’ 
reported the discovery of a possible new X-ray transient, dubbed Swift 
J023017.0+283603. The object was detected in Swift observation 
00014936012, which had taken place between 08.19 and 08.46 UTC; 
that is the notification was produced 5.2 hafter the observation. (Most 
of this latency related to the timing of ground station passes and the 
ingesting of the data by the Swift Data Center. The observation data 
were received at the UK Solar System Data Centre at 13.51 UTC.) No 
catalogued X-ray source was found at this position. Further, Swift had 
already observed this location on 11 previous occasions (the observa- 
tion target was SN 2021afkk, 4.3’ away from this serendipitous tran- 
sient), for a total of 9.6 ks. These observations had been analysed as a 
stacked image in LSXPS (Dataset ID 19690); no source was present near 
the position of SwiftJO230, with a3oupper limit of 1.5 x 10° counts s”. 
The peak count rate of Swift JO230 in the new observation was 
2.7406 x 10-? counts s” (1o errors), significantly above this upper limit, 
clearly indicating that a new transient had been discovered. 

Due tothe very soft spectrum and coincidence with the nucleus of 
the galaxy 2MASXJ02301709+2836050, this was originally interpreted 
as a TDE?”, and a high-urgency target-of-opportunity request was 
submitted to Swift. 

In the following analysis, we assume that H, = 67.36 kms? Mpc", 
Q, = 0.6847 and Q,, = 0.3153 (ref. 4). 


Observations and data analysis 

Swift follow-up observations began at 16.07 UTC on June 22 
(TO + 0.67 d). Due to Swift’s Moon observing constraint, subsequent 
observations were not available until day 4 (June 26). Daily observations 
of 1ks exposure were obtained with Swift until day 12. A subsequent 
target-of-opportunity request was submitted (Principal Investigator, 
Guolo) for weekly monitoring of this source, which began on day 21 
(2022 July 13). The initial observation showed that the source had 
resumed emitting X-rays, but in the following observations on days 27 
and 35, it was again below the detection threshold. To better quantify 
the duty cycle, we submitted regular target-of-opportunity requests 
(Principal Investigator, Evans) for daily 1 ks observations, which ran 
until 2023 March 19 (day 270) when the source entered Swift’s Sun 
observing constraint. Note that we did not obtain 1 ks per day, as each 
month proximity to the Moon prevented observations for 3-4 days. 
Moreover, due to the nature of Swift’s observing programme, our 
observations were sometimes shortened or completely superseded 
by other targets. 


Swift-XRT. XRT data were analysed using the on-demand tools of 
ref. 37, via the swiftools Python module (v.3.0.7). A 0.3-10 keV light 
curve was constructed, binned to one bin per observation; the soft 
and hard bands were set to 0.3-0.9 and 0.91-2 keV respectively. Obser- 
vations 00015231018 and 00015231019 overlapped in time, as did 
0001523143 and 0001523144. When this happens, the per-observation 


binning is unreliable, so we built light curves for each of these observa- 
tions individually, and then replaced the affected bins in the original 
light curves with those thus obtained. For each run of consecutive 
upper limit bins, we merged the limits into a single bin, using the Merge- 
LightCurveBins() function in swifttools, giving a better measurement 
of, or limit on, the quiescent flux. 

For each observation in which the light curve showed a detection 
(at the 3a level), we extracted a spectrum, fitting it with a black-body 
component absorbed by two absorbers. The first was a Tbabs model 
with N, fixed at the Galactic value of 1.12 x 107 cm” (ref. 38). The second 
was a ztbabs model with N,, free to vary and the redshift fixed at the 
value obtained from our NOT spectrum. From these fits, we obtained 
the 0.3-2 keV flux and (given the luminosity distance of 160.7 Mpc) 
the luminosity. The dependence of this luminosity on the black-body 
temperature, as reported in the main paper, is shown in Extended Data 
Fig. 2. We also obtained for each spectrum the conversion factor from 
the 0.3-10 keV count rate to the 0.3-2 keV luminosity, and so converted 
our count-rate light curve into luminosity. For the detections with too 
few photons to yield a spectral fit, the upper limits in the light curve and 
the bins created by merging (above), we used the conversion factor of 
7.54 x 10” erg per count obtained from the discovery observation. The 
resultant light curve is shown in Fig. 3, with the merged bins marked 
with crosses (and in green/cyan in the online version). To explore the 
rapid flux decay seen at the end of the first outburst (right-hand panel 
of Fig. 3), we rebinned the data to one bin per snapshot (using the 
rebinLightCurve() function), converting to luminosity using the same 
conversion factor (that obtained for the appropriate observation) 
for each bin. 

To determine the timescales of the outbursts (Table 1), we defined 
outbursts as being times when the 0.3-2 keV luminosity was above 
3 x 10* erg s” (based on visual inspection of the light curve). We built 
anew light curve, still with one bin per observation, but in which all 
bins were created as count rates with lo errors, rather than allow- 
ing upper limits (using the swifttools module with the argument 
allowUL=False passed to the light curve), and then identified each 
point that was inconsistent with L =3 x 10* erg s ‘at at least the lo 
level. The start and end times of the outbursts were then constrained 
to being between consecutive data points from this sample that were 
on opposite sides of the 3 x 10” erg s” line. The results are shown in 
Table 1. We created a Lomb-Scarle periodogram (using the astropy. 
timeseries.LombScargle Python package) to search this light curve 
for periodicity. We found possible peaks centred on 22.1 and 25.0 d, 
each with widths of ~1d. To determine their significance, we used 
a bootstrapping method in which we ‘shuffled’ the data, randomly 
redistributing the fluxes (and their errors) among the time bins, and 
then recalculated the periodogram. We did this 10,000 times. For each 
trial period in the periodogram, we identified the 99.7th percentile 
of power, that is the 30 significance threshold. The result is shown 
in Extended Data Fig. 1, along with the window function. These two 
peaks are both clearly above 3a in significance and not present inthe 
window function. 

We also investigated whether we could rule out short-period 
variations with our observing strategy, like those seen in the QPEs 
reported to date. We simulated a simplistic light curve based on the 
period of GSN 069 (and including alternating between slightly longer 
and shorter cycles as in GSN 069). For each snapshot in the real XRT 
light curve of Swift J0230, we determined the phase of the trial period 
and set the count rate either to 0.03 counts s” (on) or 0.001 counts s 
(off); fractional errors were set to typical values from our real light 
curve. We then constructed the Lomb-Scargle periodogram for this 
and repeated the bootstrap approach above. A strong signal was found 
at the nominal period, confirming that such a signal would have been 
easily recovered had it been present. Thus, we can be confident that 
there was no short-period modulation like that in GSN 069 present 
in Swift JO230. 
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Swift-UVOT. UVOT data were analysed using the uvotsource package 
in heasoft v.6.30. For the pre-outburst data, the location of Swift JO230 
was only in the field of view on five occasions. No sign of variability 
was found in these data, so we summed the images in each filter using 
uvotimsum and extracted mean magnitudes. In the initial discov- 
ery observation, UVOT gathered data in all filters, but no sign of the 
outburst was seen, any variability being swamped by the underlying 
galaxy emission. UVOT magnitudes from this observation and the 
prediscovery data are shown in Extended Data Table 1. Due to this lack 
of variability, subsequent Swift observations used the UVOT ‘filter of 
the day’, which rotates each day between the u and ultraviolet filters, 
to preserve the life of the filter wheel. No significant variability was 
seen. When visually examining the light curve, it is tempting to claim 
some variability in phase with the XRT data, but the magnitude of the 
variability is much smaller than the errors on the UVOT photometry. 
To further investigate, we rebinned the XRT light curve to one bin per 
snapshot (that is the same binning as the UVOT data, which has one 
exposure per snapshot) and disabled upper limits, forcing acount rate 
and loerror per bin. For each UVOT filter, we identified the coincident 
XRT data and then performed a Spearman rank correlation analysis 
between the XRT and UVOT fluxes. This did not account for the uncer- 
tainty on the count rates, and therefore was likely to overestimate 
the significance; however, no significant correlation was found at all, 
with P values between 0.1 and 0.9, and so more complex correlation 
mechanisms were not deemed necessary. We also attempted image 
subtraction, summing all UVOT images in the u filter (which showed 
the strongest signs of possible variability) during times of XRT detec- 
tion and non-detection, before subtracting the latter from the former. 
No evidence of an excess at the XRT position was seen. 


Nordic Optical Telescope. Spectroscopy of the host galaxy 2MASX 
JO2301709+2836050 was obtained on 2022 November 1 (day 132). 
A 2.4 ks optical spectrum was accumulated using the Alhambra Faint 
Object Spectrograph and Camera mounted onthe 2.56 m NOT located 
at La Palma, Spain. The spectrum was reduced using the spectroscopic 
data reduction pipeline PyNOT (https://github.com/jkrogager/PyNOT). 
We used a1.0 inch slit width and grism no. 4, covering the wavelength 
range -3,200-9,600 A at resolution AA/A = 360. The airmass during 
the observation was of the order of -1.1. The spectrum is shown in 
Fig. 2, and features prominent Ha, [N 11] and [O 111] emission lines, 
with acommon redshift of 0.03657 + 0.00002. A weak Hf line is also 
seen. The flux of weaker lines is often affected by the presence of stel- 
lar absorption in the continuum. To recover the pure nebular fluxes, 
we fitted the spectrum with the STARLIGHT software (http://www. 
starlight.ufsc.br/). STARLIGHT fits the stellar continuum, identifying 
the underlying stellar populations in terms of age and metallicity. 
By comparing the observed data with the output synthetic spectrum, 
the pure nebular continuum can be identified and the emission-line 
fluxes measured accurately. 

Based on this analysis, we could build the Baldwin, Phillips and 
Terlevich (BPT) diagram”. Such diagrams are widely adopted to iden- 
tify the level of nuclear activity ina galaxy. It exploits the ratio of nearby 
emission lines, minimizing the effects of extinction. The result is shown 
in Extended Data Fig. 3. 2MASX JO2301709+2836050 lies in the locus 
where a low-power AGN, alow-ionization nuclear emission-line region 
and star-forming dominated galaxies intersect. A secure classification 
of 2MASXJ0O2301709+2836050 is, therefore, not possible, but a power- 
ful AGN is clearly ruled out. 


Liverpool Telescope. The position of Swift JO230 was observed with 
the optical imager (10:0) of the 2 m LT on six different occasions 
between 2022 Jun 28 and 2022 Aug 26, using the griz filters (the first 
and last epoch were gri only). Images were processed using the default 
10:0 pipeline and downloaded from the LT archive. We co-added expo- 
sures and performed basic image subtraction versus Pan-STARRS1 


reference imaging using a custom IDL (Interactive Data Language) 
routine. Although a few subtracted image pairs show weak positive or 
negative residuals at the location of the -18 mag nuclear point source, 
there is no clear correlation in these residuals between filters or epochs, 
suggesting minimal contribution of any nuclear transient (or any AGN 
variability) to the optical flux at the sensitivity level of the LT images. 
The lack of any residual source at the location of SN 2020rht is unam- 
biguous in all images. Limiting magnitudes of the images (So) are given 
in Extended Data Table 2. 


Chandra. We requested a 3 ks director’s discretionary time observa- 
tion of Swift JO230 with Chandra (Proposal 23708869). We triggered 
this on day 93, when the Swift-XRT count rate exceeded the approved 
threshold of 0.02 counts s”, and the observations were obtained on day 
97 (observation ID 27470). Our intention was to obtain an accurate (arc- 
second or better) position of Swift J0230, to be able to say definitively 
whether it was associated with the nucleus of its host, the historical 
supernova or neither. Unfortunately, this observation occurred during 
the quiescent/faint part of the fifth outburst, and Swift JO230 was not 
detected. The 3oupper limit, converted to 0.3-2 keV luminosity assum- 
ing a90 eV black body witha Galactic absorber, is L < 8.0 x 10*' erg s7, 
consistent with the XRT measurements at the time (Fig. 3). 


Data availability 

All the Swift data are available via the Swift data archives provided 
in the USA (https://swift.gsfc.nasa.gov/archive/), UK (https://www. 
swift.ac.uk/archive/) and Italy (https://www.ssdc.asi.it/mmia/index. 
php?mission=swiftmastr). The target IDs are 00014936 and 00015231. 
Reduced Swift-XRT data for this transient are available at https://www. 
swift.ac.uk/LSXPS/transients/690. The Chandra data are publicly 
available via the Chandra data archive (https://cxc.harvard.edu/cda/), 
with sequence 704871 and observation ID 27470. The NOT data will 
be available through the NOT public interface after the expiration of 
the standard proprietary period. The reduced spectrum is available 
through the University of Leicester FigShare repository (https://doi. 
org/10.25392/leicester.data.c.6444296). The LT data will be available 
through the LT public interface after the expiration of the standard pro- 
prietary period. The photometry was included in this published article. 


Code availability 
The IDL code used for the Liverpool Telescope data analysis is available 
from D. Perley on request. 
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Extended Data Fig. 1| Period analysis of the XRT data of Swift J0230. The Lomb-Scargle periodogram of the per-snapshot binned XRT light curve is shown in red. 
The window function is in grey and the line marking 3-o significance as a function of period is in black. The two peaks above the 3-0 line and not corresponding to 
window-function peaks are centred on 22.1 d and 25.0 d. 
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Extended Data Fig. 2 | The 0.3-2 keV luminosity and blackbody temperature derived from spectral fits to the XRT observations in which Swift JO230 was 


detected. A Spearman rank test shows these to be strongly correlated (p-value: 4.5 x 10%). The errorbars reflect the 90% confidence intervals on the parameters, 
obtained using AC =2.706 in the spectral fitting. 
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Extended Data Fig. 3 | BPT (Baldwin, Phillips & Terlevich) diagram, showing left panel, the green line is the demarcation between pure star forming and 
galaxy type (HII star forming region, AGN, LINER, composite) asafunctionof composite star-forming/AGN regions, as prescribed by. The straight segments 
certain line flux ratios. The line ratios are [O iii] 5007 / HB versus [N ii] 6583 / separate proper AGN from LINERs. The SDSS galaxy catalogue object density is 
Ha (left) and [S ii] 6717,6732 / Ha (right). In both panels, the red solid lines are shown in greyscale and the position of 2MASX JO2301709+2836050 is marked by 
the theoretical models separating star-forming regions and AGN*°. Inthe the blue cross (errorbars are 1-0). 
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Extended Data Table 1| UVOT photometry from pre-discovery data and the discovery observation 


Filter 


(AB, discovery obs) 


Magnitude 
(AB mag) 


15.78 + 0.07 
16.36 + 0.07 
17.40 + 0.07 
18.33 + 0.08 
18.95 + 0.08 
18.85 + 0.07 


Magnitude 


15.82 + 0.10 
16.41 + 0.08 
17.41 + 0.09 
18.39 + 0.10 
18.88 + 0.11 
18.91 + 0.09 
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Extended Data Table 2 | Liverpool Telescope upper limits on emission at the position of Swift JO230 after subtracting the 
galaxy emission (AB magnitudes 


MJD-59752 g r a Zz 


7.20 
32.20 
25.15 
29.18 
33.21 
66.12 


>20.8 
aLI 
>20.5 
>20.6 
>22.0 
>21.5 
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